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Abstract The universe, with large-scale homogeneity, is locally inhomogeneous, cluster- 
ing into stars, galaxies and larger structures. Such property is described by the smooth- 
ness parameter a which is defined as the proportion of matter in the form of intergalactic 
medium. If we take consideration of the inhomogeneities in small scale, there should 
be modifications of the cosmological distances compared to a homogenous model. Dyer 
and Roeder developed a second-order ordinary differential equation (D-R equation) that 
describes the angular diameter distance-redshift relation for inhomogeneous cosmolog- 
ical models. Furthermore, we may obtain the D-R equation for observational H(z) data 
(OHD). The density-parameter Qm, the state of dark energy to, and the smoothness- 
parameter a are constrained by a set of OHD in a spatially flat ACDM universe as 
well as a spatially flat XCDM universe. By using of x 2 minimization method we get 
a = 0.81^ 2Q and Qm = 0.32^ ^ at lcr confidence level. If we assume a Gaussian prior 
of Q M = 026 ± 0.1, we get a = 0.93+°;g and Q M = 0.31+° *. For XCDM model, a is 
constrained to a > 0.80 but ai is weakly constrained around -1, where ai describes the 
equation of the state of the dark energy (px = wpx)- We conclude that OHD constrains 
the smoothness parameter more effectively than the data of SNe la and compact radio 
sources. 

Key words: Cosmology-dark energy-smoothness parameter 
1 INTRODUCTION 

Through recent several decades, we have entered into an accurate cosmological period and understood 
our universe more deeply. According to the cosmological principle, our universe is homogeneous and 
isotropic on large scale, but the deviations from the homogeneities are also observed. The corresponding 
researches are still exciting. 

Recently, there is mounting data from type la supernovae, cosmic microwave background (CMB) 
and the large scale structure suggest that the present universe is spatially flat and accelerated 
expanding. Combined analysis of the above cosmological observations support that an approxi- 
mately 26% of cold dark matter (CDM) and the other part 74% dominated by an unknown ex- 
otic component with negative pressure — the so-called dark energy — driving the current acceleration 
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and the Cardassian models dFreese et alJ d2002l) ) and soon proposed to explain the accelerating expan- 
sion of the universe. In the case of X-matter, the dark energy has the following property with an equation 
of state: 

Px = wp x , (1) 

where w is a constant independent of time or redshift. If a> = -1, it is reduced to the case of the 
cosmological constant (the ACDM model). 

However, these models could no t explain the ob s ervatio ns of our universe perfectly. Except for 
the cosmological constant problem dWeinberg et alJ dl989l) ). the deviation from the cosmological 
principle — the universe is homogenous and isotropic in large scale — needed to be considered. It is ob- 
vious that the matter in the universe is clustered into stars, galaxies and clusters of galaxies, rather 
than absolutely distributes uniformly and disperses everywhere in the space. It is also well known that 
the universe are grouped into superclusters, or perhaps filaments, great walls and voids in larger scale. 
Only on the scale larger than 1 Gpc does the universe appear smooth. Such problems may have effects 
on the distance redshift relation. Therefore a smoothness-paramete r a was int r oduce d to describe the 
proportion of the mean density p in the form of intergalactic matter dDver etalJd!973h ): 



Pint 



(2) 



where p mX is the mean density in the universe in the form of intergalactic matter, while p denotes the 
mean density of the whole universe, so a e [0, 1]. In the case of a — 0, it describes a universe with all 
the matter clustered into stars, galaxies and so on. while a — 1, it is a normal homogeneous universe. 
Generally, < a < 1 describes the universe partially in the form of clustered matter and the other in the 
form of intergalactic matter. 

The properties of angular diameter distance in a locally inhomogeneo u s univ e rse has been d i scusse d 



in Ref. ( Wei nberg et al 



d!989l)iZeldovichetal.l dl967l)iDashveski et al.l d!966l)iKavser et al.l dl997l) ). 
Later Dyer et al. established the Dyer-Roeder (D- R) equation to exp lain the distance-redshift relation 
in a univ erse with fractiona l intergalactic med ium ( Dyer et all ( 1973b ). as well as without intergalactic 
medium dDver et al.l dl972l) ). In the literature dSantos et al. ( 2008 )). by use of two different samples of 
SNe type la data, the Qm and a parameters are constrained by x 1 minimization fitting, which applies 
the Zeldovich-Kantowski-Dyer-Roeder (ZKDR) luminosity distance redshif t relation for a flat A CDM 
model. A ^-analysis, by using the 1 15 SNe la data of Astier et al. sample dSantos et al.l d2008l) ). con- 
strains the density-parameter to be Qm = 0.26^ ^ (2cr) while the a parameter is unlimited (all the 
values a e [0, 1] ar e allowed even at lcr). And, the analysis based on the 182 SNe la data of Riess et al. 
dRiess et al.l d2007l) ) constrains the pair of parameters to be fi M = 0.33*9'?? and a > 0.42 (2<x), which 
provides a more stringen t const raint because the sample extends to higher redshifts. 

In Ref. dSantos et al.l d2008l) ). Santos et al. has proposed to constrain a, Q.m and u> by the angular 
diameter distances of compact radio sources with XCDM model. Howerver, only are the a> and Qm 
parameters well constrained, but the a parameter is totally unconstrained in lcr level. 

As can be seen, neither SNe la nor compact radio sources data are capable of constraining the 
smoothness parameter. We could make use of other astronomical data to constrain the smoothness pa- 
rameter. It is also feasible to constrain inhomogeneous model by making use of the observational H(z) 
data (OHD ), which can be obtained by the method to estimate the differential ages of the oldest galax- 
ies. In Ref. dYi et alJ d2007l) ). Ze-Long Yi and Tong-Jie Zhang present a constraint on a flat FRW uni- 
verse w ith a matter component and a holographic dark energy component use OHD. In Ref. d Wan et alJ 
d2007l) ). Hao-Yi Wan et al. use OHD to constrain the Dvali Gabadadze Porrati (DGP) Universe. In 
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Ref. dLin et al.l (120081) 1. Hui Lin et al. successfully use OHD together with other observat ional data to 
constr ain the ACDM cosmology. The wiggling Hubble parameter H(z) are also studied dZhang et al.l 
d2008h . 

It can be conclude that OHD is a complementarity to other cosmological probes and may also 
present better constraint on the smoothness parameter. In this article, the parameters D. M , a and u> are 
constrained by totally 12 bins of OHD from Simon et al. (2005) and Ruth et al. (2008) in spatially 
flat ACDM universes as well as in the XCDM model. This paper is organized as follows: In Secf2] we 
review the basic origin of Dyer-Roeder Equation and the relationship between Hubble parameter and 
different cosmological models, which is performed in a inhomogeneous universe. In Secj3] we constrain 
the parameters Qm, cj and a from OHD. Discussions and prospects are presented in the third section. 

2 DYER-ROEDER EQUATION AND THE RELATIONSHIP BETWEEN ZKDR DISTANCE 
AND HUBBLE PARAMETER 

We consider a stellar object which emits a beam of light propagating throughout a space-time described 
by the metric tensor g^ lv . We can identify a null surface £ determined by the eikonal equation g^X^E y = 
along which the beam of light propagates. The direction of this light is the tangent vector of null surface 
i.e. the null geodesic = -S,^. The beam of light rays can be described by x? = (v,y l ), where x° = v is 
the affine parameter and y' (i = 1, 2, 3) indicates the three different directions of the propagation of the 

light. The vector field is tangent to the light ray congruence, kr = = -S, H , determines two optical 

dv 

scalars — 9 describing the convergence of the light and the shear parameter cr, 



#V cr = k^frffh v , (3) 



where nf = — - irj) is a complex vector that is orthogonal to W (k^m^ = 0). Since k M = -E,^, the 
V2 

vorticity which is connected with the light beam is zero, therefore the congruence of light is character- 
ized by th ese two optical scal ars, cr and 9. And these two optical scalars satisfy the Sachs propagation 
equations dZhang etalJdl966h l: 

6 + e 2 + \o-\ 2 = -^'R tiV kfk Y , (4) 

& + 29cr = - -C jmX nfk v m r tf , (5) 

where a dot denotes the derivative with respect to v, and % are the Ricci tensor and Ricci scalar 
respectively, and C^a is the Weyl tensor which is zero in a conformally flat FRW space-time. One 
can see that if the shear cr is initially zero, the Weyl tensor could always be zero, which automatically 
satisfies the condition in FRW space-time dPemianski et al.l d2003h ). Therefore, assuming that the light 
beam has no shear — cr — 0, we may describe the convergence and divergence of this beam of light by 
the parameter 0(empty beam approximation). 

The relative rate of the change of an infinitesimal area A on the cross section of the beam can be 
described by the optical scalar 9 (and the distortion by cr), the only parameter that characterizes the 
congruence of light, which relates to A by 

1 A 

9= . (6) 

2 A y ' 

Substi tuting Eq. d4]i into the above expression, one can reduce the optical scalar equation to dSachs et alJ 
(fl96lfo 

VA + )-<Ru V kVk v VA = 0. (7) 



The Einstein field equation is 



2 

1 $7rG 

%v - TgftvR - Agfw = 5~ r Aiv, (8) 

z c 
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Fig. 1 (color online). Angular diameter distance Da(z) as a function of redshift z for a fiat 
ACDM model. Several selected values of Qm and a are shown. Da(z) is in units of c/Hq. 



multiply each side by WW , the two terms vanish while leaving the following form, 

8nG, 



(9) 



The universe, though locally inhomogeously distributed, is homogeneous and isotropic on the largest 
scale in the mean, so we choose Robertson-Walker metric 

ds 2 = c 2 df 2 - a\t)dcr 2 , (10) 

where dcr 2 describes the spacial part of the metric and a(t) is the scale factor of the universe. We set 
a(t) — 1 at the present time and choose proper affine parameter v so that, dSchroedinger et aD dl956l) ) 

- = -^, (ID 
dr H(,a 

where ao and Ho are the present values of a and the Hubble constant, respectively. Then we have 

dr dr Hoa 

We consider a pressureless matter dominant universe in which the energy-momentum tensor has only a 
nonzero 0-0 part (with comoving coordinates), i.e. Too = P and 7^ = 0. In addition, if a/ao = (1 + z)~' 
and p/po = (1 + z) 3 , we get 

/ \2 2 

= hr 1 OP = ^+z) 5 apo- (13) 
\H a) Hi 

Substituting it into Eq. (O and then substituting the resulting equation into Eq. (0, we can finally obtain 

Va + ^q<Qm(1+z) 5 Va = 0. (14) 

Here, we use the density-parameter fi M instead of po- Due to a relationship between the angular diam- 
eter distance D& and A D k = VI" dSchneider et al.l (ll988al)jSchneider et al.l dl988bbjBartelmann et all 
dl99ll)jWatanabe et al.l dl992h ~). the Eq. flUl becomes the Dyer-Roeder Equation dDyer et al.l d 1973h >: 

r> A +^a£l M (l+z) 5 D A = 0. (15) 
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Redshift z 


H{z) 


lcr interval 


Data 


0.05 


75.4 


±2.3 


♦ 


0.09 


69 


±12 


* 


0.17 


83 


±8.3 


★ 


0.27 


70 


±14 




0.40 


87 


±17.4 




0.505 


96.9 


±6.9 


♦ 


0.88 


117 


±23.4 




0.905 


116.9 


±11.5 


♦ 


1.30 


168 


±13.4 


* 


1.43 


177 


±14.2 




1.53 


140 


±14 




1.75 


202 


±40.4 





Table 1 Observational H(z) data (OHD). The data marked with stars is from Simon et al. 
sample, and the data marked with diamonds is from Ruth et al. sample. 



where the dots denote the derivatives with respect to affine parameter v. 

It is necessary to mention that the smoothness-parameter mus t be different at various epoch of the 
univer se due to the theory of formation of the large scale structure dSantos et alJ (l2008htefstathiou et alJ 
For very high redshift, the matter in the universe must be more smoothly distributed compared 
to that of present. In this point of view, we would have to identify the smoothness-parameter a as a 
function of z, a{z) in Eq. (flSl l. especially when discussing the properties of the angular diameter distance 
at high redshift. However, because the samples of compact radio sources, SNe la and OHD are mostly 
located at low redshifts (z < 2 for Hubble parameter), we set a as a constant in the following discussion. 
We will not consider the variations of a with respect to z also because the data are neit her adequate 
nor pr ecise enough . The redshift dependence of a was discussed by Santos et al. in Ref. dSantos et alJ 
d2008h l. 

Considering Eq. (TT3T > again, we may change the variable by substituting redshift z for the affine 
parameter v and obtain, 

dz 2 d 2 D A d 2 zdD A 2 5 
(—y—- r + — ; -7- + -aCl M (l+zyD A =0. (16) 
dv dz z dv z dz 3 

Note that the univer se discussed is spatially fiat, i.e. = 0, so that Q A = 1 - fi M . Finally, after the 
substitution of variable dDemianski et al.l (l2003b ). we get a second-order ordinary differential equation 
in which the angular diameter distance D A is the function of redshift z, and D A is in the unit of c/Hq: 

d2z) A + ^ + <2 DA=0 , (17) 



dz 2 dz 



Here the initial conditions 



£>a(0) = 0, 
dD A 



dz 



1. (18) 



is satisfied. 

Besides, the functions V and Q read 



2£1m(1 + z) 3 + 3^(1 - O m )(1 + z) 3w+3 



f = 2 



^m(1+z) 4 + (1-^m)(1+z) 3w+4 

<3 



|ord M + ^(1 - Om)Ti + Z) 3 " (19) 



O m (1 +z) 2 + (1 -£2 m )(1 +z) 



3a>+2 ' 
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Fig. 2 (color online). Hubble parameter H(z) as a function of redshift z for a fiat ACDM 
model with selected values of Q.m and a. The data set in Table. [TJare also shown. 



dD A ._. 

The numerical results of Da and (hereafter DAz)) are shown in FigUJwith iterative calculations 

dz * 

by the fourth-order Runge-Kutta scheme(see Sec l3.U . From the well kn own Etherington p r incipl e — the 
relation between the luminosity distance and angular diameter distance jEtherington et alj dl933l) ). 

D L = (l+z) 2 D A , (20) 

we can get the luminosity distance as a function of z- 

At present, by the aid of the method based on the differential age of the oldest galaxies, the Hubble 
parameter can be determined as a function of redshift. It reads: 

dz 1 

df 1 +z 

which can be measured directly by the determination of dz/dt. 

Since the comoving radial distance r(z) (in units of c/H^) in fiat geometry can be expressed as 



KZ)= T^T, (22) 

Jo E(z) 

where E(z) is the expansion rate of the universe, which relates Hubble parameter to Hubble constant Ho 
in the equation: 

H(z) = H E(z), (23) 
the angular diameter distance can be written as, 

Da = (24) 
1 +z 

Differentiating Eq. d22l) with respect to redshift z and combining with Eq. I|241l . we can get the expansion 
rate of the universe expressed by Da and D' A (z) at any redshift z: 

E(z)=- . * , - ■ (25) 

(1 +z)D'(z) + D A 
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Fig. 3 (color online). Hubble parameter H(z) as a function of redshift z for a fiat XCDM 
model with selected values of oj and a. The data set in Table. Q]are also shown. 
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Fig. 4 (color online), (a) Confidence regions at 68.3%, 90.0%, 95.4%, and 99.7% levels from 
inner to outer respectively on the (Qm, oi) plane for a flat ACDM model (without considering 
the prior on Dm)- The "x" in the center of confidence regions indicates the best-fit values 
(0.28,0.97). (b) The one-dimensional probability distribution function (PDF) p for the a pa- 
rameter, (c) PDF for the £)m parameter. 
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Fig. 5 (color online), (a) Confidence regions at 68.3%, 90.0%, 95.4%, and 99.7% levels 
from inner to outer respectively on the (Dm. &) plane for a fiat ACDM model (a prior on £2m 
considered). The "x" in the center of confidence regions indicates the best fit values (0.27, 1.0). 
(b) PDF for the a parameter, (c) PDF for the Qm parameter. 




Fig. 6 (color online), (a) Confidence regions at 68.3%, 90.0%, 95.4%, and 99.7% levels from 
inner to outer respectively on the (a>, a) plane for a fiat XCDM model. The "x" in the center 
of confidence regions indicates the best fit values (-1.08,0.99). (b) PDF for the a parameter, 
(c) PDF for the a> parameter. 



3 SAMPLES AND RESULTS 
3.1 The observational data of H(z) 

In order to constrain smoothness parameter and other cosmological parameters with OHD, we need to 
integrate Eq. ( fTTI ) to obtain Da(z) and D' A (z) as a function of z (FigQ]), then from Eq. d23l and Eq. d25l ) 
to get Hubble parameter as a function of redshift z theoretically. 
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Al though it is u n able to obtain the a nalyti cal solution of Eq. dTTl l dKantowski et al.1 

(119981) jKantowski et al.l d2000t)jKantowski et al.l d200ll) ). o ne can get an approxim ate expression of 
the equation which is accurate enough to use in practical (iDemianski et al.l (120031) ). It is convenient 
for controlling the precision that we integrate the Eq. (fTTI i iteratively applying the 4th-order Runge- 
Kutta scheme and then get numerical results of D A (z) and D' A {z). Furthermore, the Hubble parameter 
H(Q. M , a), a, h\ z) can be obtained with any cosmological model at arbitrary z. 

Note that OHD, consisting of two pa rts, is obtained from two different sources. For one part, from 
Simon et al. sample dSimon et al. (120051)) , we ha ve a sample of 9 bins and z € [0, 1.75], while for the 
other, from Ruth et al. sample ( Ruth et ail d2005l) ). we choose the data separated into 3 bins. In TableQ] 
we listed all the data and errors mentioned above with their sources marked, and we plotted them in 
Figf2]and Figf3] respectively. We can compare the theretical curves of H(z) and the data in different 
models (see Sec l3.2| and Sec ]3.3t . 



3.2 Constraining a and Qm for ACDM model 

Firstly, we can find the relationship between H(z) and a, Qm- In Figf2] we plot the theoretical H(z) at 
redshift z e [0, 2] according to the last section for some typically selected a and Qm, where a flat ACDM 
model is assumed. The figure include the mean values of OHD in redshift bins and their error bars. One 
can see from Figf2]that the curve of H(z) strongly depends on Qm, i.e., the larger for Qm, the faster for 
the growth of H(z). In contrast, the smoothness-parameter a has effects on the properties of the Hubble 
parameter mainly at high redshift. For a certain Qm, theoretical H(z) with different a seem to appear 
similar at lower redshift, but they start to differ at high redshift. 
In order to constrain a and Qm we use x 1 minimization 

12 

X 2 (Hq, a, Qm) = ^ 

;=1 

where H{Hq, a, Qm; Z;) is the theoretical expecting of Hubble parameter which is determined by Eq. 
dl7) . Eq. d25l l and Eq. d23l , and i/ b s fo) is the observational values of the Hubble parameter with errors 
cr(zi) in the sample. 

In the analysis, we marginalize Hubble constant Ho by inte gration over it, and assum ed a Gaussian 
prior according to the best fitting value obtained from Ref. dBonamente et al.l d2005l) ). i.e., Hq = 
76.9^ ^km s _1 Mpc~'. On the basis of the cosmic concordance from observations, we can choose to 
add a Gaussian prior on Q M optionally, i.e., Q M = 0.26 + 0.1. We investigate the minimization both 
considering this prior(Fig[2a)) and without considering it(FigHta)). In Fig|4f a) and Fig[2a) we plot the 
regions of confidence on the Qm - a plane. The contours on the confidence of 68.3%, 90.0%, 95.4% and 
99.7% are determined by two-parameter levels 2.30, 4.61, 6.17 and 11.8, respectively. 

We also plot the one-dimensional probability distrinbution functions p (PDF) of parameters Qm and 
a. In FigHJb) and FigHJc), PDFs were plotted without considering the Gaussian prior, while in Fig|3b) 
and FigJ5J c ) we considered this prior. One can see that in 90% confidence region, 0.59 < a < 1.0 and 
0.23 < Qm < 0.41 if we consider the prior. While without any prior a 90% confidence lies in the region 
of 0.42 < a < 1.0 and 0.23 < Q M < 0.54. The "x"s indicate the model with the best fitting values 
that occur at a — 1, Qm = 0.27 and or = 0.97 , Qm = 0.28, respectively. It is clear that whether or not 
we consider the Gaussian prior, the best fitting models are nearly the same with a slightly lower than 
one, corresponding to a universe with the uniform distribution of cold dark matter, and the value of Qm 
favors other observations. 



H(H , a, Q M ; zd - H obs (zj) 
crizt) 



(26) 



3.3 Constraining a and oj for XCDM model 

In this part, we constrain a and u> for XCDM model. In this case, we set the density-parameter Qm to 
be 0.28 from the best fitting results in the Sec l3.2l In Figf3]we plot theoretical H(z) with some typically 
selected a and a>, from which one can see how these two parameters modify the theoretical curve. 
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Unfortunately, we find that the curves is not so strongly denpend on the parameters as in Sec l3.2l We 
constrain these two parameters by x~ minimization 



where H(Hq, a, co; Zi) is the theoretical value of Hubble parameter while ff bs(Zi) and cr(zi) are their 
observational values and errors, respectively. Again, marginalizing the parameter Ho, we get the regions 
of confidence on the co - a plane. Then integrating over co and a in the two-dimensional probability 
function p(a, co), we obtain PDFs of a and co. The numerical results were plotted in Fig(6] 

From Fig|6ja), we find a is only mildly constrained. The best fitting point occurs at a — 0.99 and 
co = - 1 .08 which indicates that the state of equation of XCDM is approximately that of the cosmological 
constant. All the allowable values of a in [0, 1] are permitted at the 3cr confidence level. 

4 CONCLUSIONS AND DISCUSSIONS 

In this article, we study how the inhomogeneous distribution of cold dark matter affects the Hubble 
parameter at different redshift. By using OHD from Simon et al. together with Ruth et al. in a flat 
ACDM model, the smoothness-parameter and density-parameter are constrained at different confidence 
intervals. By marginalizing the Hubble constant Ho, we found that the best fitting values are a = 0.97 
and Qm = 0.28 for ACDM model. While if we assume a Gaussian prior of Qm = 0.26 + 0.1, we found 
that the best fitting values are a = 1 and Qm = 0.27. With the XCDM model, setting Qm = 0.28, we 
get the best fitting values a = 0.99 and co = -1.08. Comparing with the constraints on the smoothness- 
parameter with samples of comp act radio sources — all the values of a from to 1 are allowed at 
68.3% statistical confidence level (ISantos etal.ld2008l) ). the constraint on luminosity distance from SNe 
la gave little better results, which is 0.42 < a < 1.0 and 0.25 < O m < 0.44 at 90% confidence level 
(ISantos et alj (120081) ). In our work, from Fig |4] and Fig|5] the empty bea m (a — 0) case is ex cluded at 
even 3cr confidence level. This result is better than the two previous works (ISantos et all (1200811?!) ) which 
constrain a parameter. We can see that OHD constrains the smoothness parameter more effectively than 
both SNe la data and the data of angular diameter of compact radio sources. 

In the meanwhile, as one can see from Figf2]and Figj3] the errors of the data is too large to constrain 
a into a relatively small interval, especially at high redshift, or to let a be a somehow accurate value. The 
statistical effects may not be neglected in this analysis. In the near future, we expect better constraint on 
a from more accurate data and/or more data at high redshift. In addition, OHD, also angular diameter 
distances of compact radio sources and luminosity distances of SNe la, may also enable us to study the 
a(z) as a function of z. 

Throughout our work, we discuss the smoothness-parameter independent of space-time, i.e., p; nt + 
Pcius is constant everywhere, and a only describes the ratio of them. In a real universe the structures 
of walls and voids could not be described by such model, therefore these complex structures, which 
might be characterized by more parameters and/or more complicated models, may need more and more 
precise data. Paul Hunt et al. have discussed the case that we located in a 200-300 Mpc void with ter ribly 
low density, which is expanding at a rate 20-30% higher than the average rate dHunt et al.l d2008l) ). In 
the future work, we should consider more complex models better describe our universe, together with 
consummate physical theories, in order to find out more accurate answer. 
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